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Abstract 

The n integrals in involution for the motion on the n-dimensional el¬ 
lipsoid under the influence of a harmonic force are explicitly found. The 
classical separation of variables is given by the inverse momentum map. In 
the quantum case the Schrodinger equation separates into one-dimensional 
equations that coincide with those obtained from the classical separation 
of variables. We show that there is a more general orthogonal parametri- 
sation of Jacobi type that depends on two arbitrary real parameters. Also 
if there is a certain relation between the spring constants and the ellipsoid 
semiaxes the motion under the influence of such a harmonic potential is 
equivalent to the free motion on the ellipsoid. 


1 Introduction 

In this paper we are concerned with the motion on a n-dimensional ellipsoid under 
the influence of a harmonic potential. The problem was first posed by Jacobi in 
the nineteenth century [|l| in the context of explicitly solvable by quadratures dif¬ 
ferential equations most of them originating in completely integrable Hamiltonian 
systems. In the last decades the problem was considered mostly by mathemati¬ 
cians obtaining new results including a description of integrals in involution and 
connection with hyperelliptic curves of genus g = n, see @- 0 - However in this 
approach no explicit separation of variables was found. Quite recently we have 
obtained a new form for the n prime integrals in involution for the free motion 
on the ellipsoid in the Jacobi parametrisation, form which is very convenient for 
proving the separation of variables and solving the Hamilton-Jacobi equation |^. 

*^Electronic mail: dita@hera.theory.nipne.ro 
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The purpose of this paper is to extend our previous results in two directions: 
to find a generalisation of the Jacobi orthogonal parametrisation and the most 
general form of the harmonic potential for which the motion is still completely 
integrable. As unexpected results we obtained a new orthogonal parametrisation 
that depends on two arbitrary real parameters generalizing the Jacobi one, and 
if the spring constants ki and ^2 along two axes of symmetry satisfy the relation 
Oi ki = 02 ^ 2 , where Oi and 02 are the squares of the corresponding semiaxes then 
the motion under the influence of such a harmonic potential is the free motion on 
the ellipsoid. In our approach the separation of variables is realized by the inverse 
of the momentum map that provides an explicit factorisation into Liouville tori. 
In the standard approach the separation of variables is a rather difficult problem, 
see e.g. where the separated variables are defined as zeros of the off diagonal 
elements of the associated Lax matrix. 

In Section 2 we define a two-dimensional family of Jacobi elliptic coordinates 
and in Section 3 we prove explicitly their orthogonality. The integrals in involu¬ 
tion for the free motion on the ellipsoid are found in Section 4 and we show that 
they are a particular case of a more general form. The most general form of the 
harmonic potential for which the motion is still integrable is given in Section 5 
and the classical separation of variables is obtained in Section 6. The separation 
of the Schrodinger equation is shown in Section 7 and the paper ends with Con¬ 
clusion. In Appendix are collected a few known properties of the Vandermonde 
determinant which are used in the paper. 


2 Generalized Jacobi’s coordinates on the 
n-dimensional ellipsoid 


Usually the n-dimensional ellipsoid is viewed as a surface immersed in a n -|- 1- 
dimensional Euclidean space, defined by the equation 


^2 2 

•M' ^ Jbr\ 

— + — + ■■■ + 
CLi 0,2 


X. 


n+l 

0,n+l 


1 


where Oi, i = 1,... ,n -1- 1 are positive numbers, Oi G R+, and we suppose that 
they are ordered such that ai > • • • > a„+i > 0. We describe the n-dimensional 
ellipsoid by intrinsic coordinates and here we slightly generalize the Jacobi’s ap¬ 
proach by defining the generalized elliptic coordinates as the solutions of the 
equation 


n+l 

E 


2=1 


bj xj 
^iibi 


1 


( 1 ) 
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with bi G R*, bi 7 ^ bj, for i ^ j, i,j = 1,... ,n + 1. We rewrite (1) in the form 


. , ^ hxj XQ{X) 

^ ai{X — bi) P{X) 


( 2 ) 


where P{X) and Q{X) are the monic polynomials 


n n+1 

■p(^)=n 

i=l i=l 

Here Wj, f = 1,..., n are the Jacobi coordinates and the factor A comes from the 
fact that A = 0 is a solution of (1) when the coordinates Xi he on the ellipsoid. 
Calculating the residues on both left and right sides of equation (2) one gets 


2 QiPi) 


1 ,... ,n + 1 


( 3 ) 


If = Oj, i = 1,... ,n + 1, (3) gives the usual Jacobi parametrization of the 
n-dimensional ellipsoid. For a* = r^, f = 1,..., n + 1, the relation (3) gives the 
orthogonal parametrization of the n-dimensional sphere of radius r and in this 
case ( 1 ) has the form 


n+l 

E 


i=l 


X-bi 


0 


Our approach is sufficiently general and gives the parametrization of any quadric 
dehned by Oj G R*. 

The parametrisation (3) depends on 2n -|- 2 parameters, half of them bi G 
R*, f = 1,... , n -|- 1 , being arbitrary real numbers and for this reason the corre¬ 
sponding Jacobi coordinates Ui could not be orthogonal. However it does exist 
an orthogonal parametrization, different from the usual one, which depends on 
n + ?) parameters. In other words there is more than one orthogonal system of 
coordinates of Jacobi type and this situation could be interpreted as a hidden 
symmetry of the problem. We found that this symmetry is two-dimensional and 
the independent parameters in equation (3) leading to orthogonal coordinates 
may be chosen as a*, f = 1 ,..., n -|- 1 , and 62 , the last two being arbitrary 
non-zero real numbers with 61 7 ^ & 2 - 

To hnd the new parametrization we deduce from (3) 


and 


Xi 


Xi 



Ui 

Ui - bi 


E 


n+l .. n 

i=l l,m=l 


Ui Ur, 


n+l 

E 


xt 


(ti/ 62) 
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Taking into account the relation (3), a careful inspection shows that the function 


n+l 

Qm ^ ^ 
i=l 


Xi 


{ui - hi){Ur 


I 7 ^ m 


is a symmetric polynomial with respect to in its n — 2 variables, having degree 
equal to n — 2 and n — 1 independent coefficients. An important property is that 
these coefficients do not depend on I and m, i.e. for all / 7 ^ m, cim dehnes a single 
function; instead of n{n — l)/2 different polynomials we have only one. Imposing 
now the condition cim = 0 for / 7 ^ m we obtain n —1 equations which can be solved 
with respect to bi. So the number of parameters of the new parametrization is 
2n + 2 — (n — 1) = n + 3. 

To see what happens let consider the case n = 2; in this case the condition 
C 12 = 0 is equivalent to the relation 


h 


0162 ~ ^2^1 + *^3(^1 ~ ^2) 

di — (I 2 


By iteration we obtain for n = 4 the additional constraint 


_ 0162 - 02^1 + 04(^1 - ^2) 

04 —- 

Ol — O 2 

Thus we suppose that the general case is given by 


0162 - 02^1 + an{bi - 62) „ , , 

bn = -, n = 3,..., n + 1 

Ol — 0,2 

and in the next Section we prove that the parametrisation (3) with bi, i = 
3,... ,n + 1, given by (4) and bi and b 2 , 7 ^ b 2 , arbitrary real non-vanishing 

numbers furnishes n orthogonal elliptic coorddinates. Before doing that we in¬ 
troduce a more uniform notation by dehning two new parameters 


a = 


Oi 62 ~ 0-2 bi 

di — d2 


and (3 


bi - &2 
d\ — d2 


Solving with respect to bi and 62 we hnd 


bi = a + P tti, i = 1,... ,n + 1 


(4') 


In the next Section we prove that the two-dimensional family of Jacobi coor¬ 
dinates whose parametrisation is given by the relation (3) with bi as in (4') is 
orthogonal for all real values of a G R and /? G R*. 
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3 Orthogonality property 


For proving the orthogonality of the new coordinates we write equation (2) 
another form 


n+l 

E 


i=l 



ai{X - bi) 


Q(A) 

P(A) 


in 


In the last relation we substitute 6 j = a + /3aj,i = l,...,n + l, and after minor 
transformation we get 


n+l 

E 


2=1 



a + (3 tti — X 


If {X-a)Q{X) \ 

/? V PW ) 


For I ^ m we have 


n+l n 

^ {ui -bi){Um-bi) 


'n+l 9 

X 


-?/ — h- E—/ 


n+l 2 

xj 


Um \ 3 U-iYi bi 

\ *=1 


i=l 


ui - bi 


( 5 ) 


Substituting = a + /? Oj in the last relation and using (5) and the property 
Q{ui) = 0 , i = 1 ,..., n, we hnd that = 0 for / 7 ^ m, i.e. the new coordinates 
are orthogonal, and our assumption expressed by the relation (4) is true. 

For I = m we have 


n+l 9 1 n+l 9 

{ui - biY dz^ (z- bi) 

and we need a calculation of the last sum. It can be obtained by differentiating 
(5) with respect to A. We use this result to write the Lagrangean in the form 


C 


1 

2 


n+l 


2=1 


n 


i=i 


( 7 ) 


where Qj = — Q'i.'^j) = dQ{x)/dx\x=uj- For a = 0 and /3 = 1 

one recovers the usual result 0 . 

Following the standard procedure we hnd the Hamiltonian of the free motion 
on the n-dimensional ellipsoid 

n n 

'H = ^PjUj-C =-2l3^g^pl ( 8 ) 

i=i i=i 


where = P{ui)/{ui — a) Q'{ui) and Uj, pj are canonical coordinates. 

Unlike the classical result we have obtained that the Hamiltonian of the 
geodesic motion is not uniquely dehned, it depends continuously and non-trivially 
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on two arbitrary real parameters a ans /3. By changing these parameters one 
changes the classical state of the system if the latter one is dehned as a point in 
the phase space but the form of the energy does not change. This property could 
be interpreted as a gauge symmetry of the classical Jacobi problem. Under the 
transformations 

x‘^i at , bi ^ a + (3 tti, a e R, f3 eR* 

the Lagrangean C and the Hamiltonian Ti remain invariant, and ipso facto the 
equations of motion. Since with given initial conditions the physical motion is 
unique the only freedom we have in proving the uniqueness is the reparametriza- 
tion of time. Generally speaking this fact rises a new problem, namely that of 
hnding all the metrics which lead to the same physical motion; in other words 
how large the hidden symmetry is, or how many geometries describe the same 
physical process. A solution to this problem could be of interest in the study of 
more complicated models arising in classical and quantized held theories. 


4 Integrals in involution 

We introduce now n integrals in involution that are linear independent. With 
this aim we dehne the symmetric functions of the polynomials Q'{uj) = Qj{uj) 


n —1 


Q 






hi) 


k=0 


( 9 ) 


where = (—1)^ ..., Wj-i, Wj+i, ■ ■ ■, Un) and cxfc is the symmetric polyno¬ 

mial of degree k. By construction the coordinate uj does not enter the symmetric 
sum \ k = 0,1,... ,n — 1. The following functions 

n 

Hk = ^ k = l,...,n (10) 

i=i 

with Hi = —H up to a numerical factor are n integrals in involution for the 
geodesic motion on the ellipsoid. 

The integrals of motion H^, k = 1,... ,n in the above form were found by us 
in ref. for another approach see e.g. Moser |^. 

An inspection of (10) shows that for each degree of freedom the contribu¬ 
tion to the Hamiltonian is given by a product of two factors. The hrst one, 
the “kinematical” factor, depends on a special structure, in our case the Van¬ 
dermonde structure dehned by the ratio /i(mi, ... ,Un) = Q'{ui), and the 

second one, the kinetic energy, f 2 {ui,pi) = pf P{ui)/{ui — a) depends on the 
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“physics”, in our case the geometry of the body. The important issue is the fac- 
torizatin fi{ui,..., Un) ■ / 2 ('Wi, Pi) where fi has no momentum dependence and /2 
depends only on a pair of canonical variables and nothing else. 

Let g{p,u) = T-[{p,u) be an arbitrary function depending on the canonical 
variables p and u which is invertible with respect to the momentum p. As we 
will see later the invertibility condition is necessary for the separation of variables 
in the Hamilton-Jacobi equation. In particular we may suppose that T-[{p,u) is 
a one-dimensional Hamiltonian. For each n G N we dehne an n-dimensional 
integrable model by giving its n integrals in involution 




n ^{i) 

5 FT)= 


( 10 ') 


We will prove the involutive property in the more general case (10'). We have 


{h,,h,} = 5 ^ 


fdH, 

dill 

dUk 


n 

- V 

dg{pj, Uj) 

( dJ 


qU) 

5hA 

\duj 

dpj 

dPj 

dujj 

i=i 

dPj 


duj 

Q'{Uj) 

duj J 


1 dg{pj,Uj) / (j) d 


9pj 




Q'{ui 


qO) ^ 11 1 

I 9(p*.“*)717; 


Q'{Ui 


1 9g{pj,Uj) d f g{j)i,Ui){S\J _^^) \ 
dpj duj Q'{ui) j 

The last step was possible because the symmetric functions and do not 
depend on Uj. Looking at the last expression it is easily seen that the partial 
derivative with respect to Uj vanishes for i = j. For i j we have to show that 


n cii) qij) qij) qO) 
duj Ui — Uj 


but this is a consequence of the following identities 


•9 Q{i) _ q{i,j) 


and 4-1 - 4-1 = («* 


u 




ihj) 
2 


where the upper index {i,j) means that the corresponding expression does not 
depend on both Ui and uj. In this way we have shown that H;} = 0 
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5 Harmonic potential 


In the following we want to find the most general form of the harmonic potential 
for which the motion on the ellipsoid under the influence of this potential is still 
integrable. Moser says that “the motion on an ellipsoid under the influence of 
a potential \x\^ is also integrable [and] this was shown already by Jacobi [^”. 
On the other hand Arnold makes a stroger statement: “Jacobi showed that the 
problem of free motion on an ellipsoid remains integrable if the point is subjected 
to the action of an elastic force whose direction passes through the center of the 
ellipsoid |@”, which might be understood as suggesting that the spring constants 
on different axes are different. We start with the most general form for the 
harmonic potential 

.. n+l 

u = -'^hxl 

i=l 

with ki 7 ^ kj for i ^ j and we look for conditions on ki for which the motion is 
integrable. We show that the most general form for U depends on the semiaxes 
of the ellipsoid and two arbitrary parameters which can be taken ki and ^ 2 - By 
substitution of the relation (3) in the above formula we get 


U 


n+l 


2 ^ 

1 = 1 




/n+l 


ki (Xi 


k=Q \ i= 


X^P'ibi 

\ 1=1 ' 




n—k 


where Sk are up to a ± sign the symmetric functions of ui,..., Un- This means 
that W is a polynomial of degree n in the elliptic coordinates dehned by the 
equations (3)-(4'). The motion described by the Hamiltonian H = 'H+U is 
completely integrable iff the coefficients of all the products Ui ■ ■ ■ ui vanish such 
that U in the new variables should have the form U = A + B{ui + ■ ■ ■ + «„)• The 
vanishing of these coefficients leads to the relations 


ki 


aia 2 {ki — /C2) ~ tti^ttiki — 112^2) 
~ Ol) 


3,4,..., n + 1 


which shows that the spring constants do not depend on the previous introduced 
parameters a and (3. Like the preceding case we dehne two new parameters 


0102(^1 — /C2) 1 j- ^2^2 ~ aiki r ^ t3 

7 =- and 0 =-, 7 , 0 G R 


®2 ~ 

such that ki has the dependence 


®2 ~ 


kii^Cli) 5 ^ 1, . . . , Tl “|- 1 


that can be viewed as the action of a special element of G (R) which transforms 
the right half-space into itself. For 7 = 0 one recovers the classical result. Thus 







we obtained that there is a two-dimensional family of coefficients /cj, which depend 
also on the semiaxes of the ellipsoid, for which the motion is integrable. Making 
all the calculation one hnds 

1 / A A ” 

W = - I n ^ + 7 + (5(ai + aa) gfc - ^ 

V ^ k =3 ^ i=l 

which is the most general form of the harmonic potential for which the motion is 
completely integrable. If 5 = 0 or ai/ci = 02/^2 the harmonic potential in elliptic 
coordinates reduces to a constant, i.e. the motion is the free geodesic motion, 
fact which is noticed for the first time. 


6 Separation of variables 


The Hamiltonian of the problem is H = 1-L + IA and we do not know yet how the 
other prime integrals look. We shall neglect the constants terms appearing in U 
such that 


H 


-2/3 Z 


PjUj) 2 

{ui - a)Q'{ui) 


5 


i=l 


(11) 


The Hamiltonians Hk, see the relations (10), depend on the symmetric functions 
A similar sum appears also in the potential, namely Si = — where 

Sk are dehned similarly to equation (9) by the relation 


Q{.X) =W{X -Ui) Sn-k (12) 

i=l k=Q 


It is easily seen that dSk/duj = For the other prime integrals we dehne 

the potentials Uk = Sk ioi k = 2,... ,n such that the integrals in involution 
are given by 


Hk 


-2/3 


i=n 

E 


P{Ui 


^ (ui - a)Q'{ui) 


2 3 

Pi + 


k = 1, 


,n 


The involution property follows straightforward 


(13) 




/ dUj dHj _ dHj dUj \ ^ ^ P{ui) pi / (q (q _ (q (q 

^ [dui dpi dpi duJ ^ (ui-a)Q'(ui) I 

The next important point is the separaton of variables. 


0 


9 



Let ~ T* (R"^) be the canonically symplectic phase space of the dynamical 
system dehned by the Hamilton fnnctions (13). We dehne the momentnm map 
by 


^ ^ R" : Mh = : Hi = -hi, i = I,... ,n}, hi eK (14) 

This application is snch that T“^(Mh) realizes the separation of variables giving 
an explicit factorisation of Lionville’s tori into one-dimensional ovals. Onr goal 
is to constrnct explicitely the application T“^(Mh) and for doing that we write 
the system (13) in a matrix form. With the notation f{u,p) = P{u)p‘^/{u — a) 
(13) is written as 


2(3 




S^V0'(«2) 

sP/Q'(u2j 




(5 


/ Si \ 

S2 

\ Sn / 


Mnltiplying to left by the matrix 


St^/Q'{u^) 


\ 


/ /(wi,Pi) \ 
f{u2,P2) 




V f{Un,Pn) J 


/ hi \ 

h2 

\k / 


/ 


.. 

• ^ \ 


n—1 

n—2 

. 1 


U2 

U2 

V 

n-1 

n-2 

iln • • 

• 1 / 


(15) 


we get a diagonal matrix with its non-zero elements eqnal to nnity that mnltiplies 
the vector colnmn (/(mi,Pi), f{u 2 ,P 2 ), ■ ■ ■,, f{un,Pn)Yi where t means transpose. 
In the proof one make nse of the relation A2 given in Appendix. The solntion is 


f{Ui,Pi) 


■—m" -t- — 
4/?2 * ^ 2(3 


k=l 


k ^ 


n—k 


i = 1,... ,n 


( 16 ) 


where the hrst term on the right hand side is a resnlt of the property Q{ui) = 0, 
see eqnation (12), written in the form 
have the classical form 


= — X]fc=o Sn-k- The above relations 


10 


ip{xi,pi, hi,... ,hn) = 0, i = 1,. 


,n 


which for hi = Ci, i = 1,... ,n give an explicit parametrization of Lionville’s tori. 
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Taking into account the form of f{ui,pi) and defining R{ui) = u^+ 


2/3 


ELi hk U. 


n—k\ 


we get 


Pi 



i = 1,... ,n 


where e* = ±1. By the substitution pi dS/dti the Hamilton-Jacobi equation 
separates and the solution has the form 


>S*(h, Ui , . . . , Un 



R{w) 

P{w) 


dw 


On the last expression one can see that all the subtleties of the problem are 
encoded by the hyperellitic curve y‘^ = P{x)R{x) whose genus is g = n. 

The motion described by the prime integrals (10') is also separable and 


n—1 

g{ui,Pi) ='^ hn-kUi, i = 

k=0 


For obtainaing the geodesic equations g{u,p) has to be invertible with respect to 
p. With the notation 7l{u) = EEo ^n-kUi the momentum is given by 

Pi = g~^{n{ui)) 

where g~^ denotes the inverse of g with respect to p and the solution of the 
Hamilton-Jacobi equation has now the form 


S(h,Ui, ...,Un 



g ^{TZ{w))dw 


The above formulae allow us to choose new canonical variables [Q, and in the 
last case we may take Qk = dik) k = 1 ,... ,n and the canonically conjugated 
variables Vi, i = 1,... ,n. The Hamilton equations are 


Qi = 0, i = l,...,n 


and therefore Qi 

gi, hi eR,i = 1 ,. 


Vi = -l,Vi = 0, i = 2,... ,n 

hi,i = 1,... ,n and Vi = -t + gi, Vk = gk, k = 2,.. 
, n. Because 


dS _ dS 
dQi dhi 


/ ’(g-'Y(n{w))w’'-‘dw 


2 = 1 


n with 
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where {g = dg ^{z)/dz one obtains the system of equations 


— 16ij + bj 


^ r'^j 

W/ (9-‘)'(K(«))i”-'A, j = 




which represents the implicit form of the geodesics, and shows that the canonical 
equations are integrable by quadratures. In the above formulae we singled out 
the hrst prime integral Tii] if we start with Tik as Hamiltonian then the change 
in the above formulae is 6ij —>• 6kj 


7 Quantisation 

For the beginning we consider the Hamiltonian Tii given by equation (13). It 
is well known that because of the ambiguities concerning the ordering of u and 
p we must use the Laplace-Beltrami operator 0. Its general form is A = 
'^{y/99^^Pj)y hj = ■ ,n, where g = det{gij)aiid gij is the metric tensor. In 

our case gij = and = P{ui)/{ui - a)Q'{ui). Let 14 denote 

the determinant of V, Eq.(15), and the determinant of the matrix obtained 
from V by removing the last row and the j-th column. Using the relations A.3 
and A.4 from Appendix we hnd that up to an inessential numerical factor 


9 = 




Ui — a 


p{ui) 

1=1 ^ ' 


and using it the Schrodinger equation generated by 7ii is written, after some 
simplihcation, in the form 




2=1 




Since the factor 14^!_^i does not depend on Ui it can be pulled out of the bracket 
and the preceding equation takes the form 


- P(n.) d / P{u,) avk 


^ ^ \l Ui - a du, 


2=1 


Ui — a dui 


4/3 




2=1 


Now we make use of the Jacobi identities A.5 and hnd 

n 


2=1 


P{u^) d 


P{ui) ddJ 


n—1 


Ui — a dui 


Ui — a dui 


4:f3‘ 




k=0 


= 0 
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which is equivalent to n equations of the form 


P{Ui) d 
Ui — a dui 





\l/i = 0, f = 1 ,... ,n 
(17) 


Here Ci = —Ei/2j3 and the other are arbitrary. 

The direct approach, starting with equation (16), is simpler the problem being 
one-dimensional and one gets the same equation (17). It has the advantage 
that the arbitrary coefficients are identified to = —hk/2(3, i.e. are the 
eigenvalues of the Hamiltonians Tik- 

In this way the solving of the Schrodinger eqnation was reduced to the solving 
of a Sturm-Lionville equation whose general form is 


d 

dx 



+ v{x)f{x) = \r{x)f{x) 


and the above equation has to be resolved on an interval [a, b]. It is well known 
that its eigenfnnctions will live in a Hilbert space iff p{x)r{x) > 0 on [a,b]. If 
p{x) has a continuous first derivative and p{x)r{x) a continnous second derivative 
then by making the following coordinate and fnnction transforms 

^ = J ^ ^ iriu)piu)ypf{u) (18) 

we bring the preceding equation to the standard form 

where 

= {p{u)r{u)YP 

p{cp) r{u) 

and u = u{(p) is the solution of the inverse Abel problem (18). 

In our case, Eq.(17), the transformation is 







and the Schrodinger equation gets 


dif'^ 


TMj, 




(19) 


where p{ip) = {R{u{(p)YP and in R{u) we made the rescaling hk hk/hi,k = 
1,... ,n, i.e. the solving of (17) is equivalent to solve (19) which represents the 
motion of a one-dimensional particle in the potential generated by R{u{(p)). 
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For n = 1 and 5 = 0 (19) is nothing else than the equation for the one- 
dimensional rotator 


dip'^ 


= 0 


with the solution ^ stc. In all the other cases we have 

to make use of the theory of hyperelliptic curves and/or 6'-functions in order to 
obtain explicit solutions. This problem will be treated elsewhere. 


8 Conclusion 

In this paper we have obtained the most general form of the harmonic potential for 
which the motion of a point on the n-dimensional ellipsoid is completely integrable 
and we found another form of the integrals in involution. The advantage of 
our approach is that separation of variables is very easy being an immediate 
consequence of the Stackel structure appearing into equations of motion. We have 
shown that the Vandermonde structure is enough powerful to allow construction 
of new n-dimensional completly integrable models. Two such models could be 
given by the one-dimensional Hamiltonians, g{u,p) = {sinu/u)p^ and g{u,p) = 
tgue^"^, see (10'). These models are interesting since in the first example g{u,p) 
is a function wich have a denumerable number of zeros and the second one has 
a denumerable number of zeros and poles, in both cases the hyperelliptic curve 
being of inhnite genus. Thus these models show that the dimension n of the 
system has no direct connection with the number of zeros and/or poles of the 
function g{u,p). 

Other examples of n-dimensional Hamiltonians are obtained for example from 


the many-body elliptic Calogero-Moser model []^ or the elliptic Ruijenaars model 
[Q, starting with the one-dimensional Hamiltonians Hcm{u,p) = p^/2 +v'^Vriu) 
and Hr{u, p) = cosh{ap) -^/l — 2{ai')‘^Vr{u) respectively, where Vriu) is the Weier- 
strass function and using the above machinery. 

In conclusion we discovered a new method for obtaining n-dimensional com¬ 
pletely integrable systems starting with one-dimensional Hamiltonians. Also in¬ 
teresting is the existence of a class of orthogonal metrics a la Jacobi for which 
the Lagrangean is gauge invariant. This result rises the problem of description 
of all the metrics that lead to the same physics. A step in this direction could be 
the revivihcation of the techniques discovered by Stackel, Levi-Civita, Painleve 
and many others which failed into undeserved oblivion. 
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A Appendix 


Herewith we collect a few known properties of the Vandermonde determinant, 
the novelty being their presentation in the context of algebraic duality. Let 
Oj, i = 1,..., n be n real (complex) numbers. We define the polynomial 

n n 

P{x) = l[{x -ai)=J2 A.l 

i=l k=0 

where Sk = (—..., a^) and ak denotes the elementary symmetric poly- 
nomilas of degree k. From A.l we have 


P{ai) =0, i = 1,... ,n 


i.e. the vectors Xi = {So, Si,..., Sn) and X 2 = (a",...,a, 1) are orthogonal 
{Xi,X 2 ) = 0 under the usual Euclidean scalar product. In order to see a few 
interesting duality relations we dehne the polynomials 


P,(x) = 


P(x) 


X 


n—1 

fc =0 


iJ) ^n-k-l 
Jb 


, J = 1, 


,n 


From the property Pj{ai) = P'{ai)5ij we deduce that the vectors Vj = {Sq \ ..., 
and Ui = ..., 1) are bi-orthogonal, i.e. 


(Vj, Ui) = P'{ai)6ij, i,j = l,...,n 


A.2 


showing that Vj and Ui are dual each other. By construction Vj does not depend 
on ttj. 

It is easily seen that the Vandermonde determinant has two dual equivalent 
dehnitions 


Vji{(ll, ■ ■ ■ , dn) 


1 

1 

1 


1 

1 

1 

Cll 

02 

dfi 

— 

s[^^ 

^ . . 

. s{^^ 


02 “^ . . 

n-1 


c(l) 

^n-1 

c(2) 

o„_i . . 

qin) 

^n-1 


A.3 


Let be the determinant obtained by removing the jth column and the 
last row in Vn, then the following relations hold 

Vni = Yl (ai-ak) 
kbj¥=l 


n Fh = (Vn) 
i=i 


A.3 
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n—l 




n 


A.4 


i.e. vjf\ is the Vandermonde determinant of the indeterminates Oi,..., a^i, bnt 


ttj. 


We give now the most general form of direct and dual Jacobi identities. By 
replacing the last row of the first form of 14 by the row , a^) and expanding 

over this row we find the identities, see e.g. II 

0 fc = 0,1,..., n — 2 

14 k = n — 1 

Ya=i ai k = n 

By replacing the last row of the dual form of 14 by , ■ • •, we hnd the dual 
Jacobi identity 




2=1 




0 /c = 0,..., n — 2 
14 k = n — 1 


A.5 


The above formula is a consequence of a more general result. Let Aij be the 
minor of the (i, j) element of the determinant 14 , i-e. 


From A.2 we get 



1 

1 

1 

1 


Oi . 

dj-i 

4+1 ■ ■ ■ 

CLn 

Aij 

0*4^ • 

.. o';; 

fl*-i 

4+1 • • • 

^2—1 


. 

.. af_\ 

fl*+i 

4+1 • • • 

n^+1 



■ 4-1 

a«-i 

4+1 ■ • • 

n-1 


n—l 

= -E 

qU)„ n-k- 

4 

S i = 1, 

.. ,n 


k=l 


and substitute it in the last row of A^j. Afterwards we multiply the first row by 
the second by S^}_ 2 i ^^e ith one etc., and add them to the last row 
and we get 


1 

1 

1 

1 

ai 

4-1 

4+1 

dn 

^i—1 

^2—1 

^2—1 

— 1 

Qj^ . . . 

4-1 

4+1 

“n 

4+^ 

4+1 

4-1 

4+1 


■ ■ ■ 


— ( 7 ® 

+ l ■ ■ ■ 

— ( 7 ® 
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= (-1)” * • • •, aj-i, aj+i, • ■ ■, ^n) 

where is the determinant obtained from Vn by deleting the last row and the 
j-column. 

For obtaining a dual result we denote by Bij the corresponding minor obtained 
from the dual form of Vn- Like the preceding case we use the relation 


q(l) 

^n-1 


n—2 


fc =0 


j = l,...,n 


and substitute it in the last row of Btj. Multiplying the hrst row by a" the 
second by etc., and adding to the last row we hnd 


1 

1 

1 

1 


^0-1) 

^(J+1) 

sin) 

c(l) 

• 

1 

c-b+i) 

qin) 



^i-1 

0 ( 1 ) 

cO-l) 

cO+l) 

q(n) 

‘^i+1 

c(l) 

-aj 

n-i-l cO'-l) 
-aj 

n-i-l qiJ+l) 

-aj 

n-i-l qin) 
-aj 


( Qj) • • • ) CLj-li • • • ) ®n) 
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